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Abstract 

We numerically calculate the exponent for the disorder averaged 
and fixed-sample decay of the energy-energy correlator in the (/-state 
random-bond Potts model. Our results are in good agreement with 
a two-loop expansion around q = 2 recently found from perturbative 
renormalisation group techniques, fulfill the correlation length bound 
v > 2/d, and give further evidence against replica symmetry breaking 
in this class of models. 

1 Introduction 

The g-state random-bond Potts model is an interesting framework for ex- 
amining how a phase transition is modified by quenched disorder coupling 
to the local energy density. For q > 2 such randomness acts as a relevant 
perturbation and for q > 4 it even changes the nature of the transition 
from first to second order (see Ref. M for a review). In the regime where 
(q — 2) is small, a score of analytical results have been obtained from the 
perturbative renormalization group, and the various expansions for the cen- 
tral charge and the multiscaling exponents for the moments of the spin-spin 
correlator compare convincingly to recent numerical work ||. 

A particularly useful way of carrying out these simulations is to consider 
the finite-size scaling of the Lyapunov spectrum of the (random) transfer 
matrix, thus generalizing the method commonly applied to the eigenvalue 
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spectrum in a pure system |3| . A definite advantage over the more traditional 
technique of Monte Carlo simulations || is that the transfer matrices allow 
for a representation in which q can be regarded as a continuously varying 
parameter |], ||, and in particular one can study small non-integer values of 
(9-2). 

The outcome of applying this method to the energetic sector of the trans- 
fer matrix, however, led to contradictory results j|. Most notably, the ex- 
ponent Xi describing the asymptotic decay of the disorder-averaged first 
moment of the two-point function (e(xi)e(x 2 )) oc \x 1 — x 2 \~ 2Xl seemed to 
be a rapidly decreasing function of q, in sharp disagreement with an exact 
bound on the correlation length exponent, v > 2/d 0, which in our notation 
reads Xi>l. 

More recent numerical work has emphasized the importance of crossover 
behavior |5j from the random fixed point to, on one side, the pure Potts model 
and, on the other, a percolation-like limit Q in which the ratio R = K\/K 2 
between strong and weak couplings tends to infinity It became clear that 
while the fixed ratio R = 2 employed in Ref. [[| seems to have been adequate 
for studying the spin sector when (q — 2) is small, in general higher values of 
R are needed to measure the true random behavior in the regime q > 4 H. 

These findings were put on a firmer ground when it was realized || that 
Zamolodchikov's c-theorem || can be used to explicitly trace out the critical 
disorder strength R*(q) as a function of q, by scanning for an extremum of 
the effective central charge. In conjunction with an improved transfer matrix 
algorithm in which the Potts model is treated through its representation as 
a loop model, this allowed the authors of Ref. || to produce very accurate 
results for the central charge and the magnetic scaling dimension in the 
regime q ^> 4. 

On the analytical side, the perturbative expansions for the first three 
moments of the energetic two-point function have been known for quite some 
time [r0|. It was however only very recently that Jeng and Ludwig |!~T 



succeded in generalizing these computations to a general iVth moment of the 
energy operator (e(xi)e(x 2 )) N oc \xi — x 2 \~ 2Xn , yielding 

X N = N(l- ^(3N -4)(q- 2f + 0(q - 2) 3 ) . (1) 

In particular this makes available the experimentally relevant exponent X' Q 
describing the typical decay of the energy-energy correlator in a fixed sample 



2 



at criticality . 

In the present publication we show that by combining the methods of 
Refs. [|, the exponents X 1 and X' Q can be quite accurately determined 
numerically for small (q — 2). In particular we find X\ > 1 in full agreement 
with the correlation length bound @ , and our results lend strong support to 
the above two-loop results of the perturbative renormalisation group. 



2 The simulations 

In order to compare our results with those of the (q — 2)-expansion, while 
on the other hand staying comfortably away from q = 2 where logarithmic 
corrections are expected Ul2| , our main series of data has q = 2.5. Iterating 
the transfer matrix for a strip of width L a large number M ^> L of times, 
we examine the probability distribution of the ratio between the two largest 
Lyapunov exponents Aq, Ai [13| in terms of the free energy gap Af(L) 



jjj ln(A /Ai). We employ the loop representation of the transfer matrix 
where each loop on the surrounding lattice is given a weight n = y/q |L4|] , 
and bond randomness is incorporated by weighing the two possible vertex 
configurations by Wi and 1/wi, where Wi is a quenched random variable that 
can assume two different values s and l/s, each one with probability 1/2 ||. 
By construction, the system is then on average at its self-dual point . The 
strength of the disorder is measured by s > 1, which is related to the ratio 
between strong and weak bonds by R = K\jKi = ln(l + s^fq)/ ln(l + ^fq/ s). 
The maximum strip width employed in the study is L max = 12. 

Following Ref . || , we start by locating the critical disorder strenght s* by 
searching for a maximum of the effective central charge. To do so, we must 
be able to determine finite-size estimates c(L, L + 2) [16j with five significant 
digits, which means that the free energy fo(L) = -^ln(A ) must be known 
with seven digit precision. These considerations fix the necessary number of 
iterations to be M = 10 s . 

Our results for c(L, L + 2) as a function of L and s are shown in Table [TJ 
For a sufficiently large system size L, these data exhibit a maximum as a 
function of s, the position of which determines a finite-size estimate s*(L), 
which converges to s* as L — > oo. From the data of Table [l|, supplemented 
by improved three-point fits [§J (not shown), we extrapolate to s*(g = 2.5) = 
2.5(1). 
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s 


c(6,8) 


c(8, 


10) 


c(10,12) 


1.0 


0.63775 


0.64; 


344 


0.65404 


1.5 


0.63774 


0.64<< 


354 




2.0 


0.63770 


0.64<< 


374 




2.1 


0.63767 


0.64<< 


377 




2.2 


0.63764 


0.64< l 


379 


0.65415 


2.3 


0.63759 


0.64< l 


380 


0.65415 


2.4 


0.63754 


0.64< l 


379 


0.65413 


2.5 


0.63748 


0.64< l 


379 


0.65411 


2.6 


0.63740 


0.64<< 


377 


0.65409 



Table 1: Effective central charge c(L, L + 2) of the q = 2.5 state model, as a 
function of disorder strength s and the strip width L. 

The fluctuations in Af(L) are examined by dividing the strip into M/m 
samples, each one of length m = 10 5 , from which the first few cumulants 
of A/(L) can be determined. As discussed in Ref. ||, the exponent X' 
is related to the finite-size scaling [[17] of the mean value (first cumulant) 



of A/(L), whereas X\ is similarly determined from the sum of the entire 
cumulant expansion. In practice, the second cumulant is roughly two orders 
of magnetude smaller than the first, and higher cumulants are expected to 
be further suppressed, even though their determination is made difficult by 
numerical instabilities. We can therefore with confidence truncate the sum 
of the cumulants after the second one. 



Resulting finite-size estimates [17[ of X' and X\ are shown in Tables |2] 
and respectively. Unlike what seemed to be the situation in the magnetic 
sector, these estimates exhibit a pronounced dependence on s. Ref. || worked 
at fixed R = 2, which for q = 2.5 would correspond to s ~ 1.7, and found 
Xq < 1 for all q > 2. We see here that the correct way to extract these 
exponents is to extrapolate the s = s* data to the L — > oo limit. With the 
help of improved two-point estimates j|] (not shown) we thus obtain 

X' = 1.02(1) X x = 1.00(1), (2) 

which verifies the bound of Ref. 0. These exponents, as well as the result 
for their difference Xq — X\ — 0.015(5), are in very good agreement with the 



4 



s 


^(6) 




*o(10) 


*o(12) 


1.0 


0.9791 


0.9513 


0.9375 


0.9293 


1.5 


1.0026 


0.9754 


0.9622 




2.0 


1.0468 


1.0191 


1.0057 




2.2 


1.0662 


1.0377 


1.0238 


1.0158 


2.3 


1.0760 


1.0470 


1.0328 


1.0246 


2.4 


1.0859 


1.0563 


1.0418 


1.0333 


2.5 


1.0958 


1.0656 


1.0506 


1.0419 


2.6 


1.1057 


1.0749 


1.0595 


1.0504 



Table 2: Effective exponent Xq(L) of the q = 2.5 state model, as a function 
of disorder strength s and the strip width L. 



s 


*i(6) 


*i(8) 


^i(10) 


*i(12) 


1.0 


0.9791 


0.9513 


0.9375 


0.9293 


1.5 


1.0016 


0.9746 


0.9614 




2.0 


1.0389 


1.0124 


0.9997 




2.2 


1.0534 


1.0269 


1.0142 


1.0065 


2.3 


1.0603 


1.0338 


1.0212 


1.0133 


2.4 


1.0671 


1.0405 


1.0279 


1.0198 


2.5 


1.0735 


1.0469 


1.0343 


1.0260 


2.6 


1.0798 


1.0531 


1.0405 


1.0320 



Table 3: Effective exponent X\(L) of the q = 2.5 state model, as a function 
of disorder strength s and the strip width L. 
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Q 




R* 


K 




Xt 










numerics 


theory 


numerics 


theory 


2.50 


2.5(1) 


3.3(2) 


1.02(1) 


1.023 


1.00(1) 


1.006 


2.75 


3.0(3) 


4.1(5) 


1.04(2) 


1.051 


1.01(1) 


1.013 


3.00 


3.5(5) 


4.7(10) 


1.06(3) 


1.09 


1.02(2) 


1.02 



Table 4: Numerical results for the critical disorder strength (s* or R*) and 
energetic scaling exponents (X' and Xi) as functions of q. The agreement 
with the two- loop expansion, Eq. ([!]), is good. 

(q — 2)-expansion; see Eq. ([I]). 

We have also performed simulations for higher values of q, where the 
discrepancy between Refs. and was even more pronounced, since s* is 
an increasing function of q. For q = 2.75 and q = 3, we had to increase the 
number of iterations to M = 10 9 in order to keep the error bars under control 
despite the higher disorder strength. In all cases we found good agreement 
with Ref. 0j and with the (q — 2)-expansion, at least in the range where 
the latter can be assumed to be valid. A summary of our results is given in 
Table §. 



3 Conclusion 

In summary, we have shown that the apparent violation of the correlation 
length bound [0] observed in Ref. ||] can be dismissed as a crossover effect 
due to the lack of tuning to the critical disorder strength. In conjunction 
with the results on degeneracy and descendents given in Ref. we would 
thus claim that the transfer matrix method ||, U can, at least in principle, 
be used to relate the entire Lyapunov spectrum to the operator content of 
the (as yet unknown) underlying conformal field theory. 

In particular, we have supplied convincing numerical validation of the 
two-loop expansion (|l|) for the energetic multiscaling exponents Our 



results also provide further evidence in favour of the replica symmetric ap- 
proach to the perturbative calculations, since the assumption of initial replica 
symmetry breaking leads to Xi = 1 + 0(q — 2) 3 [IBL which seems to be ruled 
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out by the results given in Table £|. 
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